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£H ■ Abstract 

, We present a code to compute the relic density of dark matter in the Next-to- 

' Minimal Supersymmetric Standard Model (NMSSM). Dominant corrections to the 

Higgs masses are calculated with nmhdecay as well as theoretical and collider con- 
straints. All neutralino annihilation and coannihilation processes are then computed 
| with an extended version of micrOMEGAs, taking into acount higher order corrections 

to Higgs vertices. We explore the parameter space of the NMSSM and consider in 
particular the case of a bino LSP, of a mixed bino-higgsino LSP and of a singlino LSP. 
\ As compared to the MSSM, neutralino annihilation is often more efficient as it can take 

I/"} ■ place via (additional) Higgs resonances as well as annihilation into light Higgs states. 

Models with a large singlino component can be compatible with WMAP constraints. 
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In any supersymmetric (susy) extension of the Standard Model (SM) with conserved R- 
parity, the lightest SUSY particle (LSP) constitutes a good candidate for cold dark matter. 
Recent measurements from WMAP [1] have constrained the value for the relic density of dark 
matter within 10% (.0945 < Qh 2 < .1287 at 2a). The forthcoming PLANCK experiment 
should reduce this interval down to 2%. It is therefore of the utmost importance to calculate 
the relic density as accurately as possible in any given SUSY model, in order to match this 
experimental accuracy. This has been done in the context of the Minimal Supersymmetric 
Standard Model (MSSM), with the publicly available program micrOMEGAs [2]. This code 
computes the relic density of the lightest neutralino LSP by evaluating the thermally av- 
eraged cross section for its annihilation as well as, when necessary, for its coannihilation 
with other SUSY particles. It then solves the density evolution equation numerically, without 
using the freeze-out approximation. 

The impact of the WMAP constraints on the parameter space of the MSSM has been 
widely studied [3]. In the mSUGRA inspired version of the MSSM, one needs rather special 
conditions among parameters to have a large enough annihilation cross-section and meet the 
constraints of WMAP. This is because the LSP is usually a bino and the main mechanism 
for annihilation is into fermion pairs and requires the presence of a light slepton. Given the 
direct constraints on light sleptons and on light neutralinos as well as on the light Higgs 
mass, this possibility is valid into a small corner of parameter space. Within mSUGRA 
models, to satisfy the WMAP constraint one must then appeal to specific processes such as 
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coannihilation processes, rapid annihilation through Higgs exchange or some non negligible 
higgsino component of the neutralino LSP. The latter condition is satisfied in the so-called 
focus point region of mSUGRA models. On the other hand, in the general MSSM with 
no assumption on GUT scale physic and free parameters at the electroweak (EW) scale, it 
is much easier to satisfy the upper bound on relic density from WMAP. All models with 
a mixed bino/higgsino or mixed bino/wino LSP annihilate much more efficiently and this 
even for heavy sfermions. In fact the annihilation is so efficient that the relic density often 
lies below the WMAP range, thus implying some additional dark matter component. Large 
higgsino component of the neutralino LSP means enhanced couplings to the Z boson and 
to the pseudo-scalar Higgs boson A, thus favouring annihilation channels into fermion pairs. 
Large wino component means more efficient couplings to gauge bosons favouring annihilation 
channels into W pairs. Furthermore coannihilation channels with a heavier neutralino and/or 
chargino are also enhanced for either higgsino or wino-like LSP's. In the limit where all 
sfermions are heavy the important parameters that enter the computation are those of the 
neutralino sector Mi, M2, /x, tan/3 as well as the pseudo-scalar Higgs mass tua- The trilinear 
couplings Aj can play a role for the relic density as they influence the mass of the lightest 
Higgs. 

It is well known, however, that the MSSM faces a naturalness problem - the so-called \i 
problem [4] - arising from the presence of a SUSY conserving mass term for the Higgs fields in 
the superpotential, fiH u Hd. The only natural values for the /1 parameter are either zero or the 
Planck scale. The first is experimentally excluded while the second reintroduces the hierarchy 
problem. The Next-to-Minimal Supersymmetric Standard Model (NMSSM) [5] provides an 
elegant solution to the \i problem via the introduction of a gauge singlet superfield S in 
the Higgs sector. Assuming the simplestjjossible scale invariant form of the superpotential, 
which contains the dimensionless \SH u Hd coupling, the scalar component S of the singlet 
acquires naturally a vacuum expectation value (vev) of the order of the SUSY breaking scale, 
giving rise to an effective fi = X(S) of order the EW scale. The NMSSM is the simplest 
SUSY extension of the SM in which the EW scale originates from the SUSY breaking scale 
only. Another nice feature of the NMSSM is that the fine tuning problem originating from 
the non-observation of a neutral Higgs boson at LEP, is less severe than in the MSSM [6]. A 
possible cosmological domain wall problem [7] can be avoided by introducing suitable non- 
renormalizable operators [8] that do not generate dangerously large singlet tadpole diagrams 
[9]. In addition to the MSSM fields, the NMSSM contains an extra scalar and pseudo-scalar 
neutral Higgs bosons, as well as an additional neutralino, the singlino. The phenomenology 
of the model can be markedly different from the MSSM [10]. The upper bound on the mass 
of the lightest Higgs state is larger than in the MSSM, up to 180 GeV [11]. Moreover, a very 
light Higgs boson (as light as a few GeV) is not experimentally excluded [12]. Similarly, a 
very light neutralino with a large singlino component may have escaped LEP searches [13]. 
All these properties could impact significantly the relic density computation in the NMSSM. 

Up to now there has been only a few studies of the relic density of dark matter in the 
NMSSM [14], although a detailed analysis of dark matter direct detection in this model was 
recently carried out [15]. In this paper we present a code that calculates the relic density 
of dark matter in the NMSSM. This code provides an interface between nmhdecay and 
micrOMEGAs. The FORTRAN code nmhdecay allows a precise calculation of the particle 
spectrum in the NMSSM, as well as a complete check of all the available experimental con- 
straints from LEP [16]. The parameters are then fed into a new version micrOMEGAs extended 
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to the NMSSM, which calculates all relevant cross-sections for neutralino annihilations and 
coannihilations and computes the relic density. 

The paper is organized as follows: we first describe the model, summarize the main 
features of NMHDECAY and discuss the implementation of the NMSSM into CompHEP/CalcHEP 
and micrOMEGAs. In section 3 we give results for typical case studies. Conclusions are given 
in section 4. 



2 Overview of the NMSSM and NMHDECAY 

Extending the calculation of the lightest neutralino relic density from the MSSM to the 
NMSSM, requires first to compute the neutralino as well as the Higgs boson masses and 
mixings in this model. To achieve this task, we used the FORTRAN program NMHDE- 
CAY [16], which conventions we will review now. We then discuss the implementation of 
the model into CompHEP/CalcHEP and micrOMEGAs, as well as the interface between both 
codes. We finally give a summary of the experimental and theoretical constraints taken into 
account. 

2.1 General Set Up 

In addition to the standard MSSM Yukawa couplings for quarks and leptons, the NMSSM 
superpotential contains two additional terms involving the Higgs doublet superfields, H u and 
H d , and the gauge singlet S 

W = XSR u H d + ^KS 3 . (2.1) 

As mentioned in the introduction, the superpotential in eq. (2.1) is scale invariant, and 
the EW scale appears only through the soft SUSY breaking terms in £ so ft, which in our 
conventions are given by 



-£ soft = m^ u \H u \ 2 + m^ d \H d \ 2 + ml\S\ 
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+(XA x H u H d S + ^kA k S 3 + h.c.) 

-- (M 2 A 2 A 2 + M1A1A1 + h.c.) . (2.2) 

We denote the Higgs vevs v u , v d and s such that 

o H uR + iH uI H dR + iH dI S R + iSi . 

" = — 71 — ' d = Vd + — 7f — ' s = s + — ■ ( 2 - 3 ) 

One can derive three minimization conditions for the Higgs vevs and use them to re-express 
the soft breaking Higgs masses in terms of A, k, A x , A k , v u , v d and s. It is also convenient 
to define the quantities 

v 

fi = Xs , v — ks , tan/3 = — and v 2 — v + v d . (2.4) 

Vd 
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2.2 Higgs Sector 

In the basis (H u r, HdR, Sr) one obtains the following mass matrix for the neutral scalar 
Higgs states: 

7W| 5ll = Mf sin 2 /? + n cot/3(A A + v) , 
A4| j22 = M| cos 2 /3 + fi tan/3(A A + i/) , 

/ m 2 \ 
M% 12 = \X 2 v 2 - -^J sin2/3 - fi(A x + v) , 

. . 3 A 2 t> 2 ,4 A sin2/3 , , „ . 
A^!,33 = ^ ^ +i/(4 e + 4i/), 

M 2 SA3 = Xv(2fi sin/3 - (A x + 2v) cos/3) , 

M 2 S>23 = Xv(2fi cos/3 - (A x + 2u) sin/3) . (2.5) 

After diagonalization by an orthogonal 3x3 matrix one obtains 3 neutral scalars hi 
(ordered in mass). Similarly, in the basis (H u j, H dI , Sj), the neutral pseudo-scalar mass 
matrix reads 

M 2 PjU = fi cot/3 ( A x + v) , 

Mi 22 — fi tan/3 (A\ + v) , 

M 2 P12 = fi(A x + u), 

. .9 X 2 v 2 sin2/3 . . 
M 2 P>m = — ^(A x + Au)-3A K u, 

M 2 P13 = Xv cos/3 (Ax — 2v) , 

M 2 P ^ = Xv sin/3(Ax — 2v) . (2.6) 
Eliminating the Goldstone mode by the rotation: 

cos/3 — sin/3 

sin/3 cos/3 I I G I (2.7) 







H dI 


H 








1 

the 2x2 pseudo-scalar mass matrix in the basis (A, Sj) has the following matrix elements 

, „ 2 A 2 t> 2 sin2/3 . , , 
Mp, 22 = g ( A + ^ ~ 3AreV ' 

A^'p 12 = Au(A A -2i/). (2.8) 

It can be diagonalized by an orthogonal 2x2 matrix P!-, yielding 2 neutral pseudo-scalars 
a, (ordered in mass). Finally, the charged Higgs mass is given by 

ml± = ^ {AxJrV) + M2w ' AV • (2 - 9) 
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In the MSSM limit (A — > 0, \x fixed) one obtains 2 quasi pure singlet states with masses 

m 2 s = u{A K + 4u) , m 2 P = -3A K u. (2.10) 
In the doublet sector, one can define 

™ 2 A = M' 2 RU = ^-(A x + v). (2.11) 
sm2p 

In the limit rriA 3> Mz, the masses of one scalar, one pseudo-scalar and the charged Higgs 
states are ~ rriA- The mass of the lightest scalar Higgs is bounded by 

m 2 hl < M J cos 2 2/5 + X 2 v 2 sin 2 2p . (2.12) 

In order to calculate the relic density, one needs the Feynman rules for triple Higgs 
interactions, in case the LSP annihilates into a Higgs pair through s-channel Higgs boson 
exchange. We give here the coupling of a non-singlet scalar hi to 2 quasi pure singlet pseudo- 
scalars P: 

9hiPP = v2X 2 v (Sa sin/5 + S i2 cos/3) + v2Xkv (S i2 sin/5 + Sn cos/3) . (2.13) 
The complete formulae for the triple Higgs interactions can be found in ref. [16]. 



2.3 Neutralino Sector 

In the basis (B, W, H u , Hd, S) the 5x5 neutralino mass matrix reads 

/ Mi Mz sm6w sin/3 -M z sm6 w cos/3 

M 2 -M z cos6 w sin/3 M z cosO w cos/3 

M z s'm9w sin/3 —Mz cos9\y sin/3 —fi —Xv cos/3 

—Mz sin9w cos/3 MzcosOw cos/3 — /i —Aw sin/3 

\ -Au cos/3 -Xv sin/3 2v ) 

(2.14) 

This matrix is diagonalized by a unitary matrix, N. The lightest neutralino LSP can then 
be decomposed as 

X°, = N U B + N l2 W + N 13 H U + N 14 H d + N 15 S . (2.15) 

Further, we shall talk about the bino, higgsino and singlino fractions of the LSP, which we 
define as N 2 ± , N 2 3 + N 2 A and N 2 5 respectively. 

In the limit A — > 0, the singlino is a quasi pure state with mass 

mg = 2z/. (2.16) 

For the relic density calculation, the coupling of the LSP to the Z boson is relevant. 
This coupling depends only on the higgsino components of the LSP and is proportional to 
N 2 3 — N 2 4 . The couplings of the LSP to the scalar and pseudo-scalar Higgs states will enter 
the computation of LSP annihilation through a Higgs resonance or t-channel annihilation 
into Higgs pairs. The Feynman rule for the LSP-scalar-scalar vertex reads 

9x ix1hi = g( N i2 - Nu tan9 w )(SaN 13 - S i2 N u ) 

+V2XN 15 (S a N 14 + S a Nu) + yfiSaiXNuNu - kN 2 ,) . (2.17) 
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The first term is equivalent to the XiXih coupling in the MSSM by replacing Sn = S22 = 
cos« and S12 = — S21 = sina while the last two terms are specific of the NMSSM. The 
second term is proportionnal to the the singlino component of the LSP while the last one 
is proportionnal to the singlet component of the scalar Higgs. The latter involves either 
the higgsino component of the LSP or its singlino component. When the LSP has a very 
small singlino component, then as in the MSSM, coupling to the scalars requires a mixed 
bino/higgsino LSP. Typically in models with a very light scalar Higgs, this one has an 
important singlet component and the coupling to the LSP is given by the last term in 
eq. (2.17). Similarly, the LSP coupling to a pseudo-scalar reads 

%?x?^ = 9(Ni2 - N u tan9 w )(P il N 13 - P i2 N M ) 

-V2XN 15 {P a N u + P i2 N 13 ) - V2P i3 (\N 13 N u - kN 2 5 ) . (2.18) 

2.4 Radiative Corrections 

Similarly to the MSSM, radiative corrections to the Higgs masses in the NMSSM can be 
relatively large [11]. As said earlier, in order to calculate the Higgs spectrum, we used the 
program nmhdecay [16]. Let us review now the accuracy with which radiative corrections 
are computed for the Higgs sector in this program. 

First, we assume that the Yukawa couplings and the soft terms are defined at the SUSY 
breaking scale Q = Msusy, corresponding to an average of the squark masses. Quantum 
fluctuations at higher scales are assumed to be integrated out through the standard renor- 
malization group evolution of the parameters from a fundamental scale like Mgut down to 
the scale Q. The effective Lagrangian at the scale Q can be assumed to be of the standard 
SUSY form plus soft terms. The full effective action then reads 

r cfr = Z i D^H^H, - V cS (H t ) . (2.19) 

i 

It is obtained from the effective Lagrangian at the scale Q by adding all quantum effects 
arising from fluctuations at scales <Q. (Here, Hi denotes all the Higgs fields, H u , H d and 
S.) These quantum effects can be classified according to powers of the various couplings, 
and powers of potentially large logarithms. 

Let us start with the corrections to the effective potential. It is somehow more convenient 
to classify the corrections to the (lightest) scalar Higgs mass, which is essentially the second 
derivative of the effective potential. At tree level, one can rewrite eq. (2.12) as 

m 2 hl ~ (g 2 + X 2 )v 2 (2.20) 

where g denotes the weak gauge couplings (we do not distinguish between large and small 
tan/3 here). 

The dominant one loop corrections to m\ originate from top, stop, bottom and sbottom 
loops. The corresponding corrections 8m\ to m\ are of order 

8 t m 2 hl ~ h\v 2 In (Q 2 /m 2 ) , 8^m 2 hl ~ h\v 2 , 

5 b m 2 hl ~ hiv 2 In (Q 2 /m 2 b ) , 5~ b m 2 hl ~ /i^ 2 . (2.21) 
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These contributions to the effective potential are computed exactly, without expanding in 
large logarithms or squark mass splittings. We also take into account the one loop pure weak 
(leading log) contributions of the order g A ln(Q 2 /M|). 

The dominant two loop corrections to the effective potential are of the form 

5 2 m 2 hl ~ (h 6 t + hta s )v 2 In 2 (Q 2 /m 2 ) . (2.22) 

Here, only the leading (double) logarithms are included, i.e. we neglect single logs as well 
as terms involving bottom and sbottom loops, proportional to powers of hb- 

Next, we review the contributions to the wave function normalization constants Z; in 
eq. (2.19). If evaluated for external momenta of 0(m t ) (the order of the Higgs masses), top 
and bottom quark loops yield contributions to Z{ of the form 

Z u ~ 1 + h 2 In (Q 2 /m 2 t ) , Z d ~ 1 + h 2 b In (Q 2 /m 2 ) . (2.23) 

After rescaling the Higgs fields so that their kinetic energies are canonically normalized, 
these effects generate contributions 5zm 2 hi which take the form 

Szml ~ (g 2 + X 2 )(h 2 + h 2 b ) In (Q 2 /m 2 ) . (2.24) 

Once the Higgs mass matrices are diagonalized, one may find eigenstates with masses 
niu substantially larger than m t . The Z factors are then evaluated at external momenta of 
0(mii), i.e. at the pole of the propagators. Hence for mg ^> m t the logarithms ln(Q 2 /m 2 ) in 
eq. (2.23) is replaced by ln(Q 2 / 'm 2 H ) , with coefficients depending on Higgs mixing matrices. 

2.5 Implementation into CompHEP/CalcHEP and micrOMEGAs 

The power of micrOMEGAs [2] is that given a set of parameters in the MSSM or in the NMSSM 
it first isolates the LSP before generating, for any given situation including those where 
coannihilations occur, all the necessary matrix elements of all relevant processes. Moreover 
an automatic procedure for looking for s-channel poles is incorporated into the program 
such that a more precise integration routine can be used in the event one is close to a pole. 
This way, any kind of annihilation or coannihilation that may be imagined is readily dealt 
with. This is possible thanks to the high level of automation based on the computation of 
cross sections through CompHEP/CalcHEP [17]. Moreover an important advantage is the ease 
with which a complicated model such as the NMSSM with its many fields and parameters, 
and therefore an extremely long list of defining Feynman rules, can be implemented in this 
package [18]. The gruesome task of having to code all possible vertices that appear in the 
NMSSM in order to generate the matrix elements is taken care of by yet another automatic 
procedure that only requires to define the Lagrangian in a very compact form, through 
multiplets and superfields. This step is performed by LANHEP [19], a program that generates 
the complete set of particles and vertices once given a Lagrangian [20]. The modification of 
micrOMEGAs to go from the MSSM to the NMSSM is done essentially through a modification 
of the model file through LANHEP. 

One drawback, though, is that CompHEP/CalcHEP only deals with tree-level matrix ele- 
ments while for the NMSSM, some parameters, especially in the Higgs sector (notably the 
lightest Higgs mass and coupling to bb) receive important radiative corrections. For the code 
to be of any use, these important radiative corrections need to be taken into account. They 
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are introduced in a gauge invariant and consistent way through an effective Lagrangian ap- 
proach. For the case at hand, to parametrize the radiative corrections to the Higgs masses 
and couplings, we write a general dimension four CP conserving effective scalar potential 
that involves the two Higgs doublets H u , H d and the singlet S as 

V iad = \i(H u H:) 2 /2 + \ 2 (H d H*) 2 /2 + \ 3 (H u H* u )(H d H*) 

+ X A (eH u H d )(eH* u H* d ) + X 5 ((eH u H d ) 2 + (6H* u H* d ) 2 )/2 

+ X{(H U H*)SS* + X s 2 (H d H* d )SS* + X s s (SS*) 2 /2 

+ \K(eH u H d )S 2 + (eH:H* d )S* 2 )/2 + \ s p (S* + S**). (2.25) 

With this effective potential, corrections to both the scalar/pseudo-scalar masses and the 
trilinear and the quartic Higgs vertices can be expressed in terms of the ten A's. The Higgs 
mass matrices eqs. (2.5,2.6) can be rewritten as 

M 2 S>U -> M% u + 2Ax v 2 sin 2 /3 -|s 2 cot/5 , 

A s 

M 2 Sj22 -> M 2 S 22 + 2A 2 v 2 cos 2 (3 ~y$ 2 tan/5 , 

X s 

M% 12 -> M| 12 + (A 3 + A 4 + A 5 ) v 2 sin2/3 + -j- s 2 , 

^1,33 ^ ^1,33 + 2(A^ + 4AP , 2 , 

•M-S 13 — >• A^| 13 + 2A* us sin/5 + A5 u s cos/3 , 

■^1,23 -^5,23 + 2 -^2 ws cos ^ + ^5 ws sin /3 (2.26) 

and 



M\ n - 


- M 2 P>11 


- 2A 5 u 2 cos 2 /3 


s 2 cot/3, 
2 P ' 


M% 22 - 




- 2A 5 u 2 sin 2 /? 


A s 

-ys 2 tan/3 , 


M\ 12 - 


-^P,12 


— A5 v 2 sin2/3 , 




M 2 P;iZ - 


+ -^P,33 


- A^ v 2 sin2/3 - 


- 16A; s 2 , 


M 2 Ptl3 - 


- -Mp,13 


— A5 vs cos/3 , 




M 2 P , 23 - 


- -Mp 23 


— A5 vs sin/3 . 





(2.27) 



(2.28) 



Finally, the charged Higgs mass eq. (2.9) can be rewritten as 



m 2 ± -> m 2 + - (A 4 + A 5 ) v 2 + ^ . (2.29) 

Starting from the results for the corrected Higgs masses and mixing angles provided by 
nmhdecay, we then solve for the A's. Having extracted these parameters LANHEP readily 
derives the corresponding scalar trilinear and quartic Higgs vertices. This is an extension 
of the procedure that was shown in detail in [21]. In the present implementation the full 
corrections from nmhdecay to the Higgs masses are included. Consequently the Higgs 
self-couplings are effectively corrected at the same leading order. As concerns these self- 
couplings, the results differ somehow from those one would obtain from nmhdecay. This is 
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due to the fact that, in NMHDECAY, only the one loop leading logarithms coming from third 
generation quarks/squarks loops are taken into account in the Higgs self couplings. This said 
the vertex which receives the largest correction, that is the one involving the lightest scalar 
hih^hi, does not contribute significantly to neutralino pair annihilation. On the other hand, 
as we will see next, the /^Qi^i vertex often enters processes of neutralino annhilation. For 
these vertices higher order corrections typically do not exceed a few percent when expressed 
in terms of corrected masses. 

2.6 Parameter Space Handling 

At the EW scale, the free parameters in the Higgs sector are (at tree level): A, k, m 2 Hu , 
m 2 Hd , m|, A\ and A K . Using the three minimization conditions of the Higgs potential, one 
can eliminate the soft Higgs masses in favour of M z , tan/3 and fi = As. We thus consider as 
independent parameters the following set of variables 

A, k, tan/3, fi, A\, A K . (2.30) 

The soft terms in the squark and slepton sector (which enter the radiative corrections in the 
Higgs sector) are also fixed at the EW scale. Since the gaugino soft masses M\ and M 2 enter 
the neutralino mass matrix, we will keep them as free parameters in our analysis. 

We set all these parameters in the program nmhdecay [16]. For each point in the 
parameter space, the program nmhdecay first checks the absence of Landau singularities 
for A, k, h t and below the GUT scale. For m^ olc = 175GeV, this translates into A < .75, 
n < .65, and 1.7 < tan/3 < 54. nmhdecay also checks the absence of an unphysical global 
minimum of the scalar potential with vanishing Higgs vevs. 

nmhdecay then computes scalar, pseudo-scalar and charged Higgs masses and mixings, 
taking into account one and two loop radiative corrections as mentioned in section 2.4, as 
well as chargino and neutralino masses and mixings. Finally, all available experimental 
constraints from LEP are checked: 

1) In the neutralino sector, we check that the lightest neutralino does not contribute too 
much to the invisible Z width (T(Z — > XiXi) < MeV if m^o < Mz/2, and that 
a{e + e~ -> x°iXf) < 10~ 2 pb if m^o +m^o < 209 GeV (i > 1) and a{e + e~ -> jgxj) < 10-1 P b 
if m^g +m^o < 209 GeV > 1). In the chargino sector, we check that the lightest chargino 
is not too light (m~± < 103.5 GeV). 

2) In the charged Higgs sector, we impose rrih± > 78.6 GeV. 

3) In the neutral Higgs sector, we check the constraints on the production rates (reduced 
couplings) x branching ratios versus the masses, for all of the CP-even states h and CP-odd 
states a, in all the channels studied at LEP (cf ref. [16] for details). 

For points which violate either a theoretical constraint (Landau Pole or unphysical global 
minimum) or an experimental constraint, a warning is issued by NMHDECAY. The Higgs 
masses and mixings calculated in NMHDECAY are then fed into micrOMEGAs. The input 
parameters of the NMSSM needed by micrOMEGAs are listed in table 1, while the input 
parameters of the standard model are specified in ref. [2]. We assume that the masses of 
the first two generations of sfermions are equal. micrOMEGAs recomputes the masses of 
charginos, neutralinos and sfermions at tree-level. In the present version the soft terms for 
sfermions are used and masses and mixings are computed at tree- level. It is straightforward 
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to replace this with some external program (such as nmhdecay) that calculates all masses 
for SUSY particles, as was done in the case of the MSSM. micrOMEGAs then calculates the 
LSP relic density, taking into account all possible annihilation and coannihilation diagrams. 
For triple Higgs vertices, the main corrections are included as outlined in section 2.5. For 
Higgs couplings to quarks, leading QCD corrections are taken into account using the running 
quark masses. 



name 


comment 


name 


comment 


hL 


A 


Mli, i=2,3 


masses of left sleptons 


hK 


K 


Mri, i=2,3 


masses of right sleptons 


tb 


tan/5 


Mqi, i=2,3 


masses of left squarks 


mu 


^ 


Mui, i=2,3 


masses of right u-squarks 


hLs 


A x 


Mdi, i=2,3 


masses of right d- squarks 


hKs 


A K 






MG1 


Mi 


Al 


t trilinear coupling 


MG2 


M 2 


Ab 


b trilinear coupling 


MG3 


M 3 


At 


t trilinear coupling 



Table 1: Input parameters for micrOMEGAs_nmssm. 



3 Results 

In this section, we discuss constraints on the parameter space of the NMSSM originating 
from the WMAP results on dark matter relic density. To avoid having to deal with a large 
number of parameters, we assume very heavy sfermions (mj = 1 TeV) and fix the trilinear 
sfermion mixing Aj = 1.5 TeV. Thus annihilation into fermion pairs through t-channel 
sfermion exchange and coannihilation with sfermions are suppressed. In the gaugino sector, 
we assume universality at the GUT scale, which at the EW scale corresponds to M 2 = 2M\ 
and M 3 = 3.3M 2 . The parameters of our model are thus A, k, fi, tan/5, A\, A K and M 2 . 
These parameters are free parameters at the EW scale. For the SM parameters, we assume 
a s = 0.118, mf olc = 175 GeV and m h {m b ) = 4.24 GeV. 

We concentrate on models which can differ markedly from the MSSM predictions, in 
particular models with tan/5 < 5 for which annihilation through a Higgs resonance is marginal 
in the MSSM. First we study the behaviour of Qh 2 as a function of /i and M 2 . Then, picking 
values for /i, M 2 corresponding to typical cases, we present contour plots in the A, k plane, 
where the difference between the MSSM (which we recover at small A) and the NMSSM 
appears explicitly. The parameters A\, A K are also critical as they affect the masses of the 
Higgs states and therefore the regions of parameter space where rapid annihilation through 
a Higgs resonance can take place. Finally, we present some choices of parameters for which 
the LSP is mainly singlino and the relic density still agrees with the WMAP constraints. 



10 



3.1 MSSM-like case: dependence in ji — M2 



To give an overview of the behaviour of the relic density in the /i, M 2 plane, in fig. 1 we 
consider a model with A = 0.1, k = 0.1, tan/3 = 5, Ax = 500 GeV, A K = 0. For this choice of 
parameters the singlino component of the LSP is small so that apart from the Higgs sector 
and the heavy neutralinos, the model is MSSM-like. Since the parameters are not specially 
tuned to encounter Higgs resonances one expects the predictions for the relic density to be 
rather similar to the MSSM. The LEP exclusion region arises from the limit on chargino 
pair production. For this choice of parameters the LEP limits on the Higgs sector does not 
play a role. Compatibility with WMAP is found in two different regions of the /1, M 2 plane, 
similar to the ones obtained in the MSSM at small to intermediate values of tan/3. 




200 400 600 800 



H [GeV] 

Figure 1: Contour plots for Vth 2 = 0.02, 0.0945 < tth 2 < 0.1287 (WMAP constraint) and 
tth 2 = 1 in the /1, M 2 plane for A = 0.1, k = 0.1, tan/3 = 5, A\ = 500 GeV and A K = 0. The 
hatched region is excluded by LEP constraints on charginos. 

The first region corresponds to a very thin band whith /x 3> M 2 , just at the boundary 
of the region ruled out by the LEP constraint on charginos. Although in this case the LSP 
is mostly a bino, efficient annihilation is possible via s-channel scalar Higgs (hi) exchange. 
However some fine-tuning is required to adjust the mass of the LSP to half the mass of hi. 
We will come back to this bino LSP scenario in section 3.2 

The second region is a band where /i > Mi = M 2 /2. There, the LSP is mostly bino with 
just enough higgsino component to annihilate efficiently into gauge boson pairs (WW, ZZ). 
This is essentially s-channel annihilation via scalar Higgs or Z exchange as well as t-channel 
chargino/neutralino exchange. The higgsino fraction necessary to obtain Qh 2 m 0.1 increases 
with the LSP mass, from 25% when m^o = 140 GeV to 50% when m^o = 400 GeV. When ji 
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and M 2 are large enough so that m^o > m t} annihilation into top quark pairs also contributes 
significantly. The onset of the top pair annihilation shows up as a kink on the WMAP allowed 
band. Finally, note that for M 2 < 200 GeV and fi « M 2 the relic density is above WMAP. 
This is because on the one hand the LSP bino component is large and on the other hand 
the gauge boson pair channel is not kinematically available. In section 3.4 we will study this 
mixed bino/higgsino case in more details. 

To the left of the last region, the LSP is mainly higgsino and annihilation into gauge 
boson pairs is efficient. For a higgsino LSP, coannihilation processes with charginos and 
neutralinos also contribute significantly. Therefore, the relic density is very small in this 
region. However, increasing the singlino content of the LSP may rise the relic density inside 
the WMAP allowed range. We will come back to this possibility in section 3.5. 

To the right of the WMAP allowed band, the LSP bino component is large while no 
Higgs state has the appropriate mass for a s-channel resonance (recall that we assume heavy 
sfermions, so annihilation to fermions through t-channel sfermion exchange is suppressed). 
Hence, for our choice of parameters, the relic density is large in this region of the /j,, M 2 
plane. However, for different choices of parameters, one may find areas where s-channel 
Higgs resonances bring the relic density down in this region. We will present such scenarios 
in section 3.3. 




H [GeV] H [GeV] 



Figure 2: a) Vth 2 as a function of [i for A = 0.1, k = 0.1, A\ = 500 GeV, A K = 0,M 2 = 
230 GeV, tan/5 = 5 (full) and tan/3 = 2,20 (dash), b) Relative contribution of the main 
annihilation channels for the case tan/5 = 5. Here VV includes both WW and ZZ channels 
and qq is the sum over all the quarks. 

To illustrate more precisely the main mechanisms at work for the LSP annihilation, we 
fix M 2 = 230 GeV and plot the relic density as a function of \x in fig. 2. The relative 
contributions of the most important channels are also displayed. At small values of fi, 
the LSP mass is below My/ and the main annihilation channel is into qq pairs through Z 
exchange. The bb channel is enhanced as it receives an additional contribution from h\ 
exchange. As \x increases, so does the LSP mass. Annihilation into gauge boson pairs rises 
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sharply and the relic density drops. The WW mode is typically 4 times larger than the ZZ 
mode. For larger /x, the LSP becomes less higgsino and the relic density increases (recall that 
the LSP coupling to gauge bosons depends only on its higgsino components). A subdominant 
contribution arises from the h\a\ channel. This mode is kinematically accesssible over the 
whole region probed since the scalar mass is « 118 GeV, the pseudo-scalar mass is 
m ai « 20 GeV and the LSP mass range is 75 < m^o < 109 GeV. For /i > 200 GeV, the 
Zhi annihilation channel becomes kinematically accessible. While this channel accounts for 
up to 10% of all annihilations, this contribution is not sufficient to bring the relic density 
within the WMAP range. 

In fig. 2a we also show the relic density for tan/3 = 2, 20. Changing tan/3 affects the LSP 
mass as well as its bino/higgsino fraction, thus having a large impact on the relic density. 
First, the LSP mass increases with tan/3. This explains why at small /i (where one can be 
below the WW threshold) the relic density is smaller at large tan/3. On the other hand, the 
bino component of the LSP increases with tan/3. Hence, for //>180 GeV (where the WW 
channel is kinematically open for all values of tan/?) annihilation is more efficient at low 
values of tan/?. 

3.2 Bino LSP, annihilation through h\ 




Figure 3: Region in the A, k plane for which 0.0945 < Vth 2 < 0.1287 (blue) for jj = 700 GeV, 
M 2 = 111 GeV, A x = 500 GeV, A K = 0, tan/3 = 5 and contour curve for Qh 2 = 1 (red). 
The theoretically/experimentally excluded regions are also displayed: LEP Higgs exclusion 
(grey), Landau pole (vertical lines), negative Higgs mass squared (horizontal lines) and non- 
physical global minima of the scalar potential (black). 
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If y 3> M2, the LSP is almost a pure bino (more than 99%), and the relic density is large 
(t-channel sfermion exchange is suppressed for heavy sfermions and annihilation through 
s-channel Z requires some higgsino component) unless annihilation proceeds through a s- 
channel Higgs resonance. In this case, the relic density is very sensitive to the mass difference, 
m hl - 2m Xl [22, 23]. For A = « = 0.1, tan/5 = 5, A\ = 500 GeV, A K = 0, y = 700 GeV, and 
M2 = 111 GeV (corresponding to the first WMAP allowed region in fig. 1) we have such a 
hi resonance. We first examine the dependence of Qh 2 on the specific NMSSM parameters 
A, k in fig. 3. In this plane, small values of k/X are excluded by either LEP Higgs searches, 
an unphysical global minimum of the scalar potential or a negative mass squared for the 
lightest pseudo-scalar. Large values of k give rise to a Landau pole and are also excluded. 
Most of the remaining parameter space is allowed by the WMAP constraints. 




Figure 4: a) Vth 2 vs A for y = 700 GeV,M 2 = 111 GeV, ft = 0.01,0.1,0.5, A x = 500 GeV, 
A R = and tan/5 = 5. b) Relative contribution of the main annihilation channels for the 
case k = 0.5. 

Next we fix k and study the variations of Qh 2 as a function of A in fig. 4. We also display 
the relative contributions of the most important channels. For A = 0.01, k = 0.1, one has 
m^o = 54 GeV and = 117 GeV so that — 2m^o = 9 GeV. A tiny higgsino fraction 
(<10~ 5 ) is then sufficcient to give the right amount of relic density (Qh 2 = 0.113). As A 
increases, decreases whereas m^o remains almost constant. Hence the mass difference 
decreases and the relic density drops sharply. When the LSP mass is above the light Higgs 
resonance there is no efficient mechanism for LSP annihilation and the value of the relic 
density shoots up. Note that for k = 0.5, Qh 2 first increases with A before dropping, due to 
a small increase of the hi mass with A, and so of — 2m^o. This corresponds to the small 
region at large k where the relic density is above the WMAP limit in fig. 3. 

The main annihilation channels are characteristics of the dominant decay modes of the 
h±. For small A, the main channels are X1X1 ~ * bb and rf. When A increases, the channel 
X1X1 ~^ a i a i becomes dominant, ai being a light quasi pure singlet pseudo-scalar. This is 
due to the fact that the hididi coupling increases with A, as can be seen from eq. (2.13). 
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On the other hand, the a± mass also increases with A (cf. M.p$$ in eq. (2.6) with A K = 0), 
so that for A>0.4 this channel becomes kinematically forbidden and the main channels are 
again bb and rf. 

3.3 Annihilation Through Resonances 

One of the new features of the NMSSM is its much richer scalar/pseudo-scalar Higgs sector 
as compared to the MSSM. Thus, one expects new possibilities for having Qh 2 « 0.1 cor- 
responding to annihilation of the LSP through a Higgs resonance. Starting from a point to 
the right of the WMAP allowed band in fig. 1 (where the relic density is large) we will now 
vary the Higgs sector parameters A\, A K to see whether it is possible to find resonances. A\ 
directly determines the heavy Higgs doublet mass tua, eq. (2.11), whereas A K influences the 
Higgs singlet masses, eq. (2.10). The neutralino sector does not depend on these parameters. 
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Figure 5: Vth 2 vs A x for /i = 300 GeV,M 2 = 300 GeV,K = 0.1, A = 0.1, A K = -50 GeV 
and tan/3 = 5. b) Masses of LSP, of scalars (dash) and pseudo-scalars (full) c) Relative 
contribution of the main annihilation channels. 



First, we fix A = 0.1, k = 0.1, \x = M 2 = 300 GeV and A K = —50 GeV and we vary A\. 
Results for VLh 2 , the LSP and Higgs masses, as well as the the main channels contributing 
to the LSP annihilation are displayed in fig 5. For this choice of parameters the LSP mass 
is m^o = 142 GeV with a bino fraction of 92%, a higgsino fraction of 7%, and a negligible 
singlino component. Thus the relic density is rather high, Qh 2 = 1.3 for A\ = 0. The 
main annihilation channel is into gauge boson pairs. Subdominant channels are into h\h\, 
Zh\ as well as fermion pairs. As A\ decreases, so does and for A\ « —250 GeV, the 
masses of the the heavy scalar and pseudo-scalar Higgs doublets, h 2 and a 2 , are such that 
m a2 = mh 2 = 2m^o. The relic density then drops sharply. The main decay modes of both a 2 
and h 2 being bb and rf, the dominant LSP annihilation channel near the resonance is into 
fermion pairs. When A\ w —300 GeV and one moves away from the s-channel resonance, 
the LSP still annihilates efficiently into h±ai or Zh\ through a 2 exchange. The relic density 
again falls into the WMAP allowed range. However, in this region hi is excluded by LEP. 

Next, we fix A\ = and vary A K , keeping all the other parameters as above. Results for 
Qh 2 , the LSP and Higgs masses, as well as the the main annihilation channels are displayed 
in fig 6. Note however that the A\ range in fig 6c is not the same as in fig 6a-b. The 
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Figure 6: tth 2 vs A K for // = 300 GeV,M 2 = 300 GeV,K = 0.1, A = 0.1, A x = and tan/5 = 5. 
b) Masses of LSP, of scalars (dash) and pseudo-scalars (full) c) Relative contribution of the 
main annihilation channels. 



LSP has the same characteristics as in the previous case. For A K ~ —950 GeV, the scalar 
(singlet) h 2 is such that m h2 = 2m^o, and for A R ps —90 GeV, the pseudo-scalar (singlet) 
a\ is such that m ai = 2m^o. Far from these resonances, Qh 2 is large and the preferred 
annihilation channel is into W pairs. The h 2 resonance is associated with an increase in 
the annihilation channel XiXi ~^ h\h\. However LSP annihilation through this resonance is 
not efficient enough to bring the relic density within the WMAP range. This is because the 
coupling XiXih2, eq. (2.17), is suppressed for a bino LSP and a singlet Higgs. Near the a\ 
resonance, the pseudo-scalar annihilates mainly into bb as shown in fig. 6c. Although the 
coupling, eq. (2.18), is also suppressed, annihilation is more efficient through a\ than 
through h 2 because of the p-wave suppression factor for a scalar exchange. 

3.4 Mixed bino/higgsino: fi « M 2 in the NMSSM 

We now consider in more details the mixed bino/higgsino case with /i = 220 GeV, M 2 = 
320 GeV, tan/3 = 5, A\ = 500 GeV and A K = (corresponding to the second WMAP allowed 
region in fig. 1). In fig. 7 we plot the Qh 2 dependence on the specific NMSSM parameters 
A, k. Theoretically and experimentally excluded regions are similar to those obtained in 
fig. 3. In the limit A — > 0, the singlino decouples. According to eq. (2.16), /i being fixed, the 
singlino mass is proportional to k/\. Hence, for small A and k/X > 1/3, the singlino is heavy 
and the LSP is a mixed bino (70%) and higgsino (30%), with a mass m^o = 140 GeV. Its 
relic density is within the WMAP bounds, the main annihilation channel being into WW, 
as explained in section 3.1. Smaller values of k/\ lead to a singlino LSP whose relic density 
is very large unless one has recourse to some special annihilation mechanism, as will be 
explained in section 3.5. In the singlino region, one can already see two clear resonances 
bringing the relic density down for A ~ 0.1 and k ~ 0.01. One corresponds to an h 2 resonance 
and the other to a Z resonance. However, these are located in a region ruled out by Higgs 
searches at LEP. 

At larger values of A and k, one can see a third resonance as well as a decrease of the 
relic density below the WMAP allowed region. To see more precisely what occurs in this 
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Figure 7: Contour plots for Qh 2 = 0.02 (dash), Vth 2 = 1 (full) and region for which 0.0945 < 
Vth 2 < 0.1287 (blue) in the A, k plane for // = 220 GeV, M 2 = 320 GeV, A x = 500 GeV, 
A K = and tan/3 = 5. Theoretically and experimentally excluded regions are labeled as in 
fig. 3. 
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Figure 8: a) Vth 2 vs A for y, = 220 GeV, M 2 = 320 GeV, k = 0.2, A x = 500 GeV, A K = GeV, 
tan/? = 5 (full) and tan/5 = 2,20 (dashed), b) Components of the LSP c) Masses of LSP, 
of scalars (dash) and of pseudo-scalars (full) d) Relative contribution of main channels for 
neutralino annihilation. 
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area of the parameter space, we fix k = 0.2 and plot Qh 2 as a function of A in fig. 8. We also 
plot the LSP components, the LSP and Higgs spectrum, and the main annihilation channels. 
At small values of A, one can see from fig. 8b that the singlino almost decouples and one 
recovers the MSSM. The value of Qh 2 for A — > and tan/3 = 5, fig. 8a, is equal to the MSSM 
prediction for equivalent parameters. The main annihilation channel, fig. 8d, is WW through 
Z exchange, with subdominant Zh\ and h\h\ contributions coming from hi exchange. As A 
increases, the relic density decreases until one encounters the h% resonance for A ~ 0.3, fig 8c. 
The relic density then drops and the main annihilation channels reflect the hi decay modes, 
with a dominant d\d\ final state. Above the resonance, Qh 2 keeps decreasing as the LSP 
becomes dominantly singlino. Recall that for \x and k fixed, the singlino diagonal mass term, 
eq. (2.14) goes like 1/A while the singlino mixings are proportional to A. This result for the 
relic density is counter-intuitive, since for a singlino LSP, annihilation should be suppressed 
and relic density large. The main annihilation channel in this case is h\d\ through t-channel 
Xi exchange. The decrease of Qh 2 with A can then be understood by an increase of the 
XiXi^i and XiXi a i couplings as given in eqs. (2.17,2.18) for a singlino LSP. However, this 
area of the parameter space is excluded by Higgs searches at LEP, as shown in fig. 7. We 
will investigate in section 3.5 whether it is possible to find regions in the parameter space 
with a singlino LSP allowed both by LEP and WMAP constraints. 

A similar dependence of Qh 2 on A is observed for different values of tan/3, in fig. 8a. The 
mass of the LSP varies with tan/3, so the value of A for which resonant annihilation occurs 
is shifted. As we have discussed in section 3.1, when the main annihilation channel is into 
gauge boson pairs, annihilation is more efficient at low values of tan/3 where the LSP is less 
bino-like. For tan/3 = 2 however, no decrease of Qh 2 occurs at large values of A. This is due 
to the fact that both the XiXi a i an d XiXiai couplings are much weaker for low values of 
tan/3. 

3.5 Singlino LSP 

We explore now scenarios satisfying both LEP and WMAP constraints with a predominantly 
singlino LSP. For this we scanned over the whole parameter space of the NMSSM in the 
range A < 0.75, « < 0.65, 2 < tan/3 < 10, 100 < fi < 500 GeV, 100 < M 2 < 1000 GeV, 
< A x < 1000 GeV and < — A K < 500 GeV. We found three classes of models: a mixed 
singlino/ higgsino LSP that annihilates mainly into h\di and VV, an almost pure singlino 
that annihilates through a Z or Higgs resonance and a singlino where dominant channels are 
coannihilation ones. In table 2 we show a selection of benchmark points along these lines. 

The first scenario is one for which \i <C M 2 and the LSP is a mixed higgsino/singlino. 
We give two different examples in table 2. For point 1, the LSP is 88% singlino and 12% 
higgsino, with a mass of 122 GeV. The main annihilation mode is h\d\ through t-channel Xi 
exchange, h\ and ai being both mainly singlet (88% and 99% respectively). As we have seen 
in section 3.4 this is due to the enhanced couplings X?X? a i(^i) f° r large values of A. The 
annihilation into h\d\ being too efficient for large values of tan/3, we chose tan/3 = 2 for this 
point. Annihilation of the higgsino component into W pairs accounts for the subdominant 
channel. For point 2, the LSP mass is 148 GeV and the main annihilation channel is WW, 
which necessitates some sizeable higgsino component (29% here). The subdominant channel 
h\d\ implies the singlino component of the LSP (69%). The singlino fraction cannot be 
maximal here since there are no important annihilation channels for a pure singlino when A 
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117 


116 
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0.04 
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0.04 


m h2 [GcV] 
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114 
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0.96 


0.01 
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0.00 


0.96 


m Ql [GeV] 


114 
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p'2 
* 12 


0.99 


1.00 


0.99 


1.00 


0.99 


1.00 


Q.K 2 


0.1092 


0.1179 


0.1155 


0.1068 


0.1124 


0.1023 


Channels 


ha (73%) 
VV (13%) 
Zh (8%) 
hh (3%) 
qq (2%) 
11 (1%) 


VV (75%) 
ha (17%) 
hh (5%) 
Zh (2%) 


qq (65%) 
11 (35%) 


qq (93%) 
11 (7%) 


aa (92%) 

qq (7%) 
ll (i%) 


^ X (81%) 
xSx§ -> X (15%) 
Xlxt - X (2%) 

qq (2%) 



Table 2: Benchmark points with a singlino LSP satisfying both LEP and WMAP constraints 
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is not so large and the channel h\a\ is not dominant. 

We give in table 2 the characteristics of three typical scenarios with light singlinos. The 
only efficient annihilation mode for a very light singlino (below 50 GeV) is via a Z or a Higgs 
resonance. For point 3, the LSP is 87% singlino with a mass of 35 GeV and annihilates 
through a Z exchange. In this scenario, although the singlino component dominates, a 
12% higgsino component is sufficient to ensure efficient annihilation through the Z . As 
for the main annihilation channels they are characteristic of the Z decay modes, mainly 
into quark pairs and lepton pairs. Neutralinos of 30 GeV that satisfy the WMAP upper 
bound were also found in SUGRA models with non-universal gaugino masses [24]. Because 
the scalar/pseudo-scalar Higgs states in the NMSSM can be much below the Z mass while 
passing all the LEP constraints, we expect to find even lighter singlino LSP's in models 
where 2m^o pa m^. We show two examples (point 4 and 5) of such models in table 2. In 
both cases the LSP annihilates via a light scalar dominantly singlet (Sf 3 ~ 1). This scalar 
decays either into bb (point 4) or, when kinematically accessible, into a±ai (point 5), the a\ 
being also mainly singlet. The Higgs sector of such models is of course severely constrained 
by LEP, in particular the limit on the SM-like scalar, here the second scalar, h 2 . For this 
reason most scenarios with light singlinos have tan/5 ~ 3 which is the value for which the 
lightest visible (ie non singlet) Higgs mass, m^, is maximized [11]. Note that a light singlino 
requires k C A and not too large value for /i. The singlino masses for point 4 and 5 are 
respectively 10 GeV and 15 GeV. We found points in the parameter space with singlinos as 
light as a few GeV. 

For k < A <C 1, the LSP is heavy with a large singlino component. No efficient anni- 
hilation mechanism is then available. However coannihilation with heavier neutralinos and 
charginos can be very efficient especially for a higgsino-like NLSP. Point 6 in table 2 gives 
an example of such a scenario. The LSP is 96% singlino with a mass of 204 GeV. The mass 
difference with the NLSP x\ * s H GeV. The coannihilation channels are overwhelmingly 
dominant. The X2 higgsino component is just enough (28%) for efficient annihilation. The 
main channels are X2X2{Xi) ~^ ti,bb and correspond to annihilation through h% and a 2 ex- 
change. For this point, /13 and a 2 belong to the heavy Higgs doublet with rriA ~ 470 GeV, 
so that we are close to a (double) resonance. Such a resonance is not necessary though, in 
order to have efficient X2 annihilation. We also found points in the parameter space with a 
heavy singlino where the dominant channel was X^X^iXi) ~^ W through Z exchange. 

4 Conclusions 

In the NMSSM, the same mechanisms as for the MSSM are at work for neutralino anni- 
hilation: into fermion pairs through s-channel exchange of a Z or Higgs, into gauge boson 
pairs through either Z/h s-channel exchange or t-channel exchange of heavier neutralinos or 
charginos. The new feature of the NMSSM is the presence of additional Higgs states, which 
means additional regions of parameter space where rapid annihilation through a s-channel 
resonance can take place. We found that annihilation through a Higgs resonance is dom- 
inant in large regions of the parameter space and this even at low to intermediate values 
of tan/3. Furthermore, new channels also open up since light Higgs state can be present. 
For example annihilation channels into Zh, hh, ha or even aa can contribute significantly to 
the relic density. These proceed through s-channel Z/h/ a exchange or t-channel neutralino 
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exchange. Despite the additional annihilation channel available because of a richer Higgs 
sector, annihilation of neutralinos is not always favoured in the NMSSM. In general the 
singlino component of the LSP tends to reduce the annihilation cross-section. Therefore one 
expects that the relic density of dark matter strongly constrain models with a large singlino 
component. We found however regions of the parameter space where a singlino LSP gives 
the right amount of dark matter, either for large A, s-channel resonances into a Z or a Higgs, 
or coannihilation with Xi- 

For a bino LSP, compatibility with the WMAP result can be recovered in the NMSSM. 
This however requires tuning the parameters of the model such that 2m^o is only a few GeV 
below the mass of a Higgs boson. In models where M2 = 1M\ as we have discussed here, 
this would also mean that one has a chargino mostly wino with a mass m~± pa 2m^o pa . 
Such a degeneracy between the chargino mass and the Higgs mass could be measured at 
the ILC If the Higgs boson is mostly doublet, it might be hard to disentangle the NMSSM 
from the MSSM. If it is singlet, mass relation between the NMSSM Higgs states would be 
different from what expected in the MSSM. A mostly singlet Higgs state could eventually 
be observable at the LHC with high luminosity [25]. 

In the case of a mixed bino/higgsino LSP, annihilation relies, as in the MSSM, on s- 
channel Z exchange and t-channel neutralino/chargino exchange. If the singlino state is not 
heavy and decouples, ie A not too small and k not too large, the five neutralino states might 
be visible at the LHC/ILC. This would be a clear signature of the NMSSM. 

Finally, in the singlino LSP cannot be too large. One therefore would expect 

visible higgsinos at the LHC. The singlino LSP would however appear at the end of the 
decay chain in any sparticle pair production process, which might complicate the detection 
task as it was the case at LEP [13]. 

As a concluding remark, one should mention that, if one assumes minimal flavour struc- 
ture, b — > S7 might impose strong constraints on the NMSSM parameter space, especially 
for large values of A where the charged Higgs mass, eq. (2.9), could be lighter than in the 
MSSM. A complete study of the constraints coming from flavour physics should be the next 
step of the phenomenological study of the NMSSM. 
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